
import libnum

a1= 1285367317452089980789441829580397855321901891350429414413655782431779727560841427444135440068248152908241981758331600586
a2= 1109691832903289208389283296592510864729403914873734836011311325874120780079555500202475594

q=libnum.gcd(a1,a2)
print(q)

kq=958026509070750981060635343194230402253811613396673763765882
p=1158310153629932205401500375817
q=827089796345539312201480770649


x=a1//kq
print(x)
print(x//p+1==p)

kq=q*(p+1)

q*(p+1)*(p*(p-1)+1)
q*(p+1)*(p2-p+1)
q*(p3-p2+p+p2-p+1)
qp3+q




print((x//p+1)**3*q+q==a1)
print(kq//q)
